ASYMPTOTICS OF THE AVERAGE HEIGHT OF 2-WATERMELONS 

WITH A WALL. 



MARKUS FULMEK 

Abstract. We generalize the classical work of de Bruijn, Knutli and Rice (giving 
the asyniptotics of the average height of Dyck paths of length n) to the case of p- 
waternielons with a wall (i.e., to a certain family of p nonintersecting Dyck paths; 
simple Dyck paths being the special casep =1.) We work out this asymptotics for the 
case p = 2 only, since the computations involved are already quite complicated (but 
might be of some interest in their own right). 



1. Introduction 

The model of vicious walkers was originally introduced by Fisher jH] and received 
much interest, since it leads to challenging enumerative questions. Here, we consider 
special configurations of vicious walkers called p-watermelons with a wall. 

Briefiy stated, a p-watermelon of length n is a family Pi, . . . Pp oi p nonintersecting 
lattice paths in Z^, where 

• Pi starts at (0, 2i — 2) and ends at {2n, 2i — 2), for i = 1, . . . ,p, 

• all the steps are directed north-east or south-east, i.e., lead from lattice point 
(z,j) to (z + l,j + l) or to (z + l,j-l), 

• no two paths Pi, Pj have a point in common (this is the meaning of "noninter- 
secting"). 

The height of a p-watermelon is the y-coordinate of the highest lattice point contained 
in any of its paths (since the paths are nonintersecting, it suffices to consider the lattice 
points contained in the highest path Pp] see Figure ^ for an illustration.) 

A p-watermelon of length n with a wall has the additional property that none of the 
paths ever goes below the line y = (since the paths are nonintersecting, it suffices to 
impose this condition on the lowest path Pi; see Figure |21 for an illustration.). 

In j2], Bonichon and Mosbah considered (amongst other things) the average height 
H{n,p) of p-watermelons of length n with a wall, 

H{n,p) 



^(all p- watermelons of length n) 

X h ■ #(all p-watermelons of length n and height h), 

h 

and derived by computer experiments the following conjectural asymptotics [3 4.1]: 



H{n,p) ~ V(1.67p - 0.06) 2n + o{Vn) . (1) 
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Figure 1 . A 6-watermelon of length 46 and height 20 




Figure 2. A 3-watermelon with a wall of length 11 and height 12 




The purpose of this paper is to work out the exact asymptotics for the simple special 
case p = 2. This will be done by imitating the classical reasoning of de Bruijn, Knuth 
and Rice [S] for the case p = 1 (i.e., for the average height of Dyck paths). However, 
even the case p = 2 involves rather complicated computations. In particular, we shall 
need informations about residues and evaluations of a double Dirichlet series, which 
we shall (partly) obtain by imitating Riemann's representation of the zeta function [HI 
section 1.12, (16)]. 

1.1. Notational conventions. For k,n & "Z, we shall use the notation introduced in 
|12j for the rising and falling factorial powers, i.e. 

(n)'' := if /c < 0, {nf := 1, {n)'' := n ■ {n + 1) ■ ■ ■ {n + k - 1) if A; > 

and 

(n)- := if A; < 0, {nf := 1, (n)- := n ■ {n - 1) ■ ■ ■ {n - k + 1) if A; > 0. 
For the binomial coefficient we adopt the convention 

h)-=l^ ifO<A;<n, 
\kj ' |o else. 
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Moreover, we shall use Iverson's notation: 
[some assertion] = 



1 if "some assertion" is true, 
else. 



1.2. Organization of the material presented. This paper is organized as follows: 

• In Section El we present exact enumeration formulas for the average height of 
p-watermelons with a wall in terms of certain determinants. Moreover, we make 
these formulas more explicit (in terms of sums of binomial coefficients) for the 
simple cases p = 1 and p = 2. 

• In Section El we first review the classical reasoning for the asymptotics of the 
average height of 1-watermelons with a wall, which was given by de Bruijn, 
Knuth and Rice [S]. Then we show how this reasoning can be modified for the 
case of 2-watermelons with a wall. 

• In appendix El we summarize background information on 

— Stirling's approximation, 

— certain residues and values of the gamma and zeta function, 

— a certain double Dirichlet series and Jacobi's theta function 
which are needed in our presentation. 

1.3. Acknowledgements. I am very grateful to Professor Kratzel for pointing out to 
me how the poles and residues of certain Dirichlet series can be obtained in a simple way 
by using the reciprocity law for Jacobi's theta function, and to Christian Krattenthaler 
for many helpful discussions. 



2. Exact enumeration 
For a start, we gather some exact enumeration results. 

2.1. The number of p— watermelons with a wall. We have the following general- 
ization of the enumeration of 1-watermelons with a wall (i.e., Dyck paths) of length n, 
which is given by the Catalan numbers 

^ 1 /2n\ 

C n = 

^ ^ n + l\ n J 

Proposition 1. The number C{n,p) of all p-watermelons with a wall of length n is 
given as 

P-I /2n+2j\ 
j=0 \ n } 

Proof. This is a special case of Theorem 6 in ^2] • D 

2.2. 1— watermelons with a wall and height restrictions. In order to obtain the 
average height, we count p-watermelons with a wall of length n which do not exceed 
height h. 

To this end, we employ the following formula (see p. 6, Theorem 2]): 
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Theorem 1. Let u, d he nonnegative integers, and let b, t he positive integers, such 
that —h<u — d<t. The numher of lattice paths from (0, 0) to {u + d,u — d) , which do 
not touch neither line y = —b nor line y = t, equals 

^{{u-k{b + t)) ~ {u-k{b + t) + b)y 



Corollary 1. Let i, j and h he integers such that < 2i,2j < h. The numher of all 
lattice paths from (0, 2i) to {2n, 2j) which lie between the lines y = and y = h > is 

m(n, 1, A) := g ((„ _ , ^ ^ - („ ^^^^ + 2) + l) ) ' 

The special case i = j = can he written as 



m{n, h) := m{n, 0, 0, h) 



n + l\ n 



2n 



^\\n-k{h + 2)-l) '^\n-k{h + 2V^\n-k{h + 2) + v)' 



k>l 



Proof. Set u = n — i + j , d = n + i — j , b = 2i + 1 and t = h + l — 2iin^. □ 

Corollary 2. For n > 0, the numher of all p-watermelons with a wall of length 2n, 
which do not exceed height h, is given hy the following determinant: 

p-i 

C{n,p,h) = m{n,ij,h) . (6) 

t,j=0 

Proof. For h > 2p — 2 this follows by a direct application of the Lindstrom-Gessel- 
Viennot method [TTj . 

For < h < 2p — 2, the determinant equals (as it should), since m{n,i, j,2i) = 
and m(n, i,j, 2i + 1) + m{n, i + 1, j, 2z + 1) = for all j. □ 

2.3. The average height of 1— watermelons with a wall. The following is a con- 
densed version of the reasoning given in [S] . Denote by m(n, h) the number of all Dyck 
paths, starting at (0, 0) and ending at (2?t,, 0), which reach at least height h. By (jSl), we 
obtain 

m(n, h) = m{n, n) — m{n, h — 1) = 1) — C(n, l,h — 1) 

y-v / / 2n \ _ / 2n \ f 2n 

k>l 



n-k{h + l) -IJ \n-k{h + l)J \n-k{h + l) + l 



So, the average height of 1-watermelons with a wall (i.e., Dyck paths) of length n is 



1 " if " ' 

^'^ h=l ^'^ \ h=0 



4 + (n + l) I 2^d(A;) ), (7) 

,A:>1 Vn/ 
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where d{k) denotes the number of positive divisors of k. Introducing the notation 

( ^" ) 

S{n,a):=J2m^'^, (8) 

fc>l \n) 

we arrive at 

H{n, 1) = (n + 1) (^S{n, 1) - 2 S{n, 0) + S{n, -1)) - 1, (9) 
which is equivalent to equation (23) in jS]. 

2.4. The average height of p— watermelons with a wall. In generalization of the 
above notation, denote by m(n,p, h) the number of all p-watermelons of length n, which 
reach at least height h, i.e., 

rn{n,p, h) = C{n,p) — C{n,p, h — 1) . 

Clearly, we have the following exact formula for the average height of p-watermelons of 
length n: 

^ n+2p~2 

H{n,p) = — r V m{n,p,h). (10) 

C{n,p) ^ 

2.4.1. The average height of 2 -watermelons with a wall. Let 

f 2n W 2n \ 

S{r^, a,b) = J2i: d(g^d(j, k)) ^""^•^;;;r'^' • (11) 

i>l fc>l [n) 

From (jni), straightforward (but rather tedious) computations lead to the following for- 
mula: 



H{n.2) = -^^—L_i^^n + iyS,{n) + S,{n)J - 1, (12) 
where 

Si(n) = -20(n - l)(n + 2)S{n, 0) + 15n{n + l)(5(n, -1) + S{n, 1)) + 
(n + 3)(6^(n, -1) - 165(n, 0) + QS{n, 1)) + (n - 2){6S{n, -1) + 
8S{n, 0) + 6S{n, 1)) - 6n{n + 3){S{n, -2) + S{n, 2)) + 

(n + 2)(n + 3)(5(n,-3) + 5(n,3)), (13) 
S2(n) = ^(n, -2, -2) - 5(n, -1, -3) - 2S{n, -1, -2) + S{n, -1, -1) + 2^(n, -1, 0)- 
S{n, -1, 3) + 2S{n, 0, -3) - 4^(n, 0, 0) + 25(n, 0, 3) - S{n, 1, -3)- 
2S{n, 1, -2) + 2S{n, 1, -1) + 2S{n, 1,0) + 5(n, 1, 1) - S{n, 1, 3) + 
2S{n, 2, -2) - 2S{n, 2, -1) - 25(n, 2, 1) + S{n, 2, 2). (14) 

3. Asymptotic enumeration 

3.1. Asymptotics of the average height of 1— watermelons. In the case of 1 
watermelons, the asymptotic of the average height (0) is well-known, see pUj Proposi- 
tion 7.7] or equation (34)]: 

H{n, 1) ~ ~ ^ + *^ . (15) 
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We repeat the classical reasoning of de Bruijn, Knuth and Rice jS], in order to make 

clear the basic idea, which we shall also employ for the case p = 2 later. 

Proof of fll5|) ; From Q it is clear that we need to investigate the asymptotic behaviour 

of S{n,a) = Y^^=i^i^) (2n\'' ■ Note that the sums S{n,a) in (jH)) are multiplied with 

a factor of order 1. So if we are interested in the asymptotics of H[n,l) up to some 
0(n~"), we need the asymptotics for S{n, a) up to 0(n~°~^); for our case, a = 1 — e is 
sufficient. 

The basis of the following considerations is the asymptotic expansion of the quotient 
of binomial coefficients (see appendix IA.1|) . 

3.1.1. The asymptotics of S{n,a) for a fixed, n ^ oo. First we observe that 

( ^" ) 

\n+a-k) r^l I 2e\\ 

=0(exp(-n )) 

if l^^l ^ n^~^ ■, i.e., if /c > n^"'"'^ + a (see (I29|l in appendix lA.lj) . Therefore, the sum of 
all terms with k > n'^^'^ + a is negligible in (jH}, being 0{n~^) for all m > 0, and we 
may take ^ = o(n'-^^ in (OTHl . 

Next, we take up to order and substitute x Pulling out the leading 

_fc2 

term and expanding the rest with respect to k gives 



(5a + 4a^) k'^ k'^ ak^ 



0(n-2+^). (16) 



3 6 n-^ 3 . 

Now we consider the following function 

g{n,h):=Y,kU{k)e-''l\ (17) 

k>l 

and observe that here the terms for k > n^^^^'^ are again negligible: 

k^ d{k) e"^'/" = O (n^™) for all m > 0. 

Hence we directly obtain from (fTB|) : 
S{n,a) = 1 1 - - 1 9in,0) + i- - gin, 1) + [^^] 9{n,2) 



,' Aa^ + 5a\ , 



(This is equation (27) in (Hj.) Note that the coefficients for g{n, k) are odd functions of 
a for odd k and obtain: 

S{n, 1) - 2 S{n, 0) + S{n, -1) = -- g{n, 0) + ^ g{n, 2) + O {n-^+'g{n, 0)) . (19) 



So we reduced our problem to that of obtaining an asymptotic expansion for g{n,b). 
Note that we need this information only for b even. It follows from the computations 
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presented in appendix IA.2.11 that g{n,2b) = 0(n^~^^^'^log{n)^, and that we have for all 
m > 0: 

g{n,0) = ^v^log(n) + (^7 - ^log(2) j V^+^ + ©(n^™) , 



n /I 3 1 \ 

g{n, 2) = -V^\og{n) + ( - + -7 - - log(2) j + ©(n"'") . (20) 

Inserting this information in (fT^ we immediately obtain the desired result (fT3j) . □ 



3.2. Asymptotics of the average height of 2— watermelons. We shall modify the 
reasoning from section appropriately. In doing so, it turns out that we have to deal 
with the double Dirichlet series {k'^+py integers a, 6 > 0. Proposition El (see 

section IA.3j) states that this series is convergent in the half-plane '^{z) > a + b + 1 
and defines a meromorphic function Z{a,b] z) : C C which has a simple pole at 
z = a + b + 1, and an additional simple pole atz = a + 6+ | only if a = or 6 = 0. 
Hence, we can write 

Z(a, b; z) = """'^ 1 + c,„, + o(z-a-b-^^. (21) 

z - a - b - ^ \ I J 

Given this "implicit" definition of the numbers ^ we will show: 



E{n, 2) ~ 0m(^-2co,o + 8ci,o - 9ci,i - 9c2,o + 15c2,i 

+ 35c2,2 + 5c3,o - 35c3,i) - 2 + C>(n-i+^) . (22) 

Using the representations of the constants Ca,b by certain integrals (see ()49j) in appen- 
dix ^3)), we obtain the following approximative asymptotics by numerical integration 
(carried out with Mathematica) 

E{n, 2) ~ 2.57758 - 2 + O (^r2~5+^ j , 

which conforms well to Bonichon's and Mosbah's conjecture (Q) for the case p = 2, 
which yields approximately 2.56125 ^Jn. Figure El shows the quotient q{ri) = ^ 57^53 2 
for small n. For example, g(lOOO) = 1.00734. 

Proof of (j22|l .- Note that in (|T2|) . the "single sums" S{n, a) are multiplied with a rational 
function in n of order at most 3, while the "double sums" S{n, a, b) are multiplied with 
a factor of order 4: So if we are interested in the asymptotics of H{n, 2) up to some 
0(n~"), we need the asymptotics for S{n,a) up to 0{n~°'~^) and for S{n,a,b) up to 
0(n~"~^); for our case, a = 1 — e is sufficient. 



3.2.1. The asymptotics of S{n,a) for a fixed, n 00. Basically, we repeat the compu- 
tations from section 13.1.11 The only difference is that we need higher orders now. After 
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Figure 3. Illustration of the quotient q{n) = 2 57^8 '^-2 ' - ^^00. 
q(n) 

1 . 15 




200 400 600 800 1000 



some calculations, we obtain: 

24 (4n2 + 20n + 89) ^(n, 0) 4(96^2 + 1065^ + 3656)^(^,2) 
Si n = + 



^288^2 + 4060n + 12213) g{n, 4) ^ 8 {Sn"^ + 107n + 335) g{n, 6) 



3n' 3n* 

Recall that in appendix IA.2.11 it is proved that g{n,2b) = 0(n^"^^/2 log(n)). Moreover, 
the arguments in app endix I A . 2 . 1 1 show that we have (in addition to (HOI)) for all m > 0: 

3^2 /I 9 3 \ 

9{n,A) = — v^log(ri) + ( 2 16^ ~ g log(2) j n^V^ + 0{n-"') , 



15n 



3 / 22 45 15 \ 

-v^log(n) + ( + ^7 - — log(2) j n=^v^ + 0(n-'") 



105r2'' , — , , , /II 315 105, , , 1 . , — 
^/(^,8) = -^v^log(n) + ( Y + -^7 - — log(2) ) n^V^ + 0{n ™) 



(7(n,10) 



945n 
128 



^ / N ^689 2835 945, 5 ^ 

-V7mlog(n) + ^ 7 log(2) n-'Jim + Oin ™) . 

^ ^ V 64 128 ' 64 7 ^ ^ > 



Inserting this information in (j23j) we immediately obtain the first part of the desired 
result: 

\2 



3.2.2. The asymptotics of S{n,a,b) for a, b fixed, n — > 00. Basically, we mimic the 
computations from section 13.1.11 In doing so, we are led to consider the following 
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function 

k>l l>l 

Observe again that the terms for fc, / > n}/"^^^ are neghgible in this sum. For obtaining 
the following formula, we made use of the fact that g{n,2a,2b) = 0(n"'~^^^^^ (which is 
shown in appendix IA.2.2|1 : 

^ , , 1 / 192 1632 8736\ , , f7Q8 8928 61744\ , 

1 / 1152 18432 161488\ , , / 576 10336 99336\ , , , 
+ 2 ( + —) 2- 2> + ( — + — - — ) 4. 0) 



/384 10784 138128\ , , , 1 /512 7872 58368\ , 
+ - — + — j 2) + 2 " ~ ^ ~ j " 

128 4192 300624\ , ^ / 256 6848 1517888\ , ^ 

—^\g(n,6,0) +{-— + — _-^)^(n,6,2) 



/2432 62368\ , ^ ^, 1 / 256 416 \ , ^ ^, /544 225488\ , „ 

/398912 960\ oo^^/^31736 256 \ ^ , 1328^(r2,8,6) , 64^(r2,8,8) 

8672 64 ^ / 128 47504\ 7856^7(r2, 10, 4) 32^?(n,10,6) 

^(n, 10, 0)+ — — - -—rr gin, 10, 2) 



5^11 3^10; ' ' ' \?>n^^ 15^12; ' ' 45^13 Sn^^ 

456^(n,12,0) 2576^(n, 12, 2) 64^(n,12,4) 16^(n,14,0) 

5ni2 + 15ni3 + 9^^^^ + 9^I^3 

32o(n, 14,2) , , , 

- 9^14 ^ +Q(n-^+^^(n,0,0)). (25) 

From the results obtained in appendix lA. 2.21 and IA.3.1'1 we easily derive the following 
asymptotic expansions (using the "implicit" definition of the numbers given in (j?H) ): 

^(^,0,0) = — + ^ (^2co,o - log(n) - ^ J - 27 J + O (n"'") , 

2a, 0) = + (4c„o - log(n) - (a + - 27) ) + O (n"-) 
2-2a-2fe-3 a+f,('2a)!('2M! /7r3^ ^ Cab(2a + 26)!\ ^/ , , 

for all m > 0. Inserting the information from (j^Uj) in shows that all the log(?T,)- 
terms cancel, as well as all evaluations of the digamma function t/; (see appendix lA. 2. 1|) . 
So we obtain the second part of the desired result: 



2 

(n + l)^S2(n) 



12 (2n + 1) 

11 

/7m( — - - 2co,o + 8ci,o - 9ci,i - 9c2 



6 

- 15c2,i + 35c2,2 + 5c3,o - 35c3,i) + 0(n-i/2+.) ^ ^27) 
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Inserting the expressions (j2D) and (P7|) in (fT^ gives the desired result ((221) • D 

Appendix A. Background information and relevant results 

A.l. Stirling's approximation applied to quotients of binomial coefficients. 

We have the following asymptotic series for log(r(2;)), valid for |arg z\ < 11—6, < 6 < tt, 
\z\ 00 (see in equation (3.10.7)]): 

iog,r(.)) . (. - i) iog(.) - + !^ H- g p^yfri) ■ P«) 



where Bj denotes the j-th Bernoulli number. 
Setting X = we thus obtain 

2n 



tXy 00 00 00 \ ( 00 00 00 00 



^"+"7^^ = exp I -2n (— + — + — + — + ... ) + (— + — + — + 
[^^) ^ y V 2 12 30 56 y V 2 4 6 

(X^ + X^ + X^ + X^ + ...)+— \ \ \ h... 



6n ^ ■■■^^3^30 12 45 4 

1 /x2 x^ x^ llx« \ 1 /x2 x^ llx^ 143x* 



rv> 9 3 14 ■ ■ ■ 7 V30 4 10 40 

'5x^ 5x'^ 91x^ 65x^ 



n9 V 66 6 1^ 



+ ... \^0{fn^^) . (29) 



Note that ^'^n*' is zero for |x| > 1. The approximation given by ()29|) is very good if, 
say, |x| < \. 

A. 2. Integral representations of the exponential function and applications. 
A. 2.1. The asymptotics of g{n,b) for b fixed, n 00. Starting with the formula 

e"^ = — : / r(^) x-'dz for c> 0, x > 1, (30) 
27ri Jc-ioo 

(see [TJ (2.4.1)]) and using 

azf = j2d{k)k-^, (31) 



fc>i 



we obtain 



k>l -^^-''^ 

— / n^r(^)C(2z-6)'d^, 



as s 



2txi 

where c > Denote the integrand in the above formula by Gi{b] z) 

For any fixed positive number q and 3?(s) > — g, we have (^(s) = O j 

00. Since '^{z) becomes small on vertical lines, we can shift the line of integration to 
the left as far as we want to, if we take into account the residues of our integrand Gi{b\ z). 
There is a double pole at z = and possibly some simple poles at 2; = 0, —1, —2, .... 
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For obtaining the residues, we use the power series expansion 

^ oo 



n=l 



where 7 is Euler's constant and 7„ = hnim^oo {Ylh=i^ ^ (logO" — {n + 1) ^ (log/)"'*'^) 



(see jHl 1.12, (17)]), which gives the Laurent expansion at ^ 

C{2z - hf = ^ + , , + ■ ■ ■ . 

Combining this with the series expansions 

= n^o (1 + log(n) {z - zq) + . . .) , 
T{z) = T{zo){l + ij{zo) (z-zo) + ...), 

for zq = where ipi^z) is the digamma function (i.e., the logarithmic derivative of the 
gamma function see [01 1-7]), we can easily express the residue of our integrand 
Gi{b;z) at z= ^: 

For our purposes, we need tp{z) at positive integral or half-integral values z, which can 
be derived from the following information: 

ip{l) = —7 (see jHl section 1.7, equation (4)]), 

V'Q) = -7-21og(2) (see m P- 104]), 

n— 1 ^ 

%jj{z + n) = ip{z) + : (see [HI section 1.7, equation (10)]). 



z + ] 



The residue of Giip] z) aX z = —m is 

where Bk denotes the fc-th Bernoulli number (see [HI 1.12, (20)]). Note that this number 
is non-zero only if h is odd or m = 6 = 0. 

The sum of fl32j) and (jHS)) for all m > gives an asymptotic series for g{n, h). 

A. 2. 2. The asymptotics of g{n,a,b) for a, b fixed, n — > 00. In the same manner as in 
section IA. 2. 11 we obtain 

g(n,a,b) = y^^^^^ n^V{z)k^f(e + f)-'<lz 

1 /"C— ioo 

= — / T{z) y d(gcd(fc, /)) A;" f (k^ + 1^)'" dz, 

2vri ic-ioo 



k,l>l 



where c > ^±|±i. 
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For k and / fixed, set j = gcd{k,l). Then we may write k = kij and I = hj with 
gcd(fci,/i) = 1. This leads to 

1 ioo 

gin, a,b) = — V{z) V d(j) ^-+""2^ kll\ [kl + llY' dz 

t/C— loo 



' c— lOO „■ 1. ; 1 

gcd(fci,/i)=l 



-| /"C— ioo 

ZTTl — ' 

C— lOO _ } 



ki,h>l 
gcd(fci,ii)=l 

Now we get rid of the constraint gcd(A;i, li) = 1: 
Proposition 2. We have the following identity: 

gcd(fc,/)=l 

Proof. By inclusion-exclusion, the left-hand side equals the sum over all pairs {k, I) 
minus the sum over all pairs {k,l), where some prime number p divides gcd(/c,/), plus 
the sum over all pairs {k,l), where the product of two different primes piP2, divides 
gcd(/c, /), and so on: 

A;"/^ ^ k^^l^ ^ ^ {kpYilpf 



Etv L ^ — V rv I ^ — V — V 

(1.2 _|_ ]2\z ~ 2^ (h.2 I 12\z ^ 2-^ 2-^ 



< + /2). ^ n.2 + 12Y ^ ^ (n, \2 + n yy 

k,l>l ^ ' k,l>l ^ ' p prime fc,«>l ^' ' \ ' > 

gcd(fc,0=l 



y-v ^ {kpiP2y{lpiP2Y 



Pl7^P2 prime fe,(>l 

k^^l^ 



p prime Pi^P2 prime / fc,/>] 



(F + py 



\p prime 



^ ^ / fe,/>l ^ ^ 



The product in the last line is the reciprocal of the Euler product for C{2z — a — b) ([HI 
p. 225], which proves the assertion. □ 

Thus, we arrive at 

1 pc-ioo r,ajb 

g{n, a,b) = — I V{z) C{2z - a - b) T^-py (35) 

t/C— 100 J V ^ / 



k,l> 



Denote the integrand in the above formula by 6*2(0., b; z). 

Again, we may shift the line of integration to the left as far as we want to, if we take 
into account the residues of our integrand G2(b; z). Computing the poles and residues 
clearly depends on some information about the double Dirichlet series involved. This 
information will be provided in the next (and final) subsection. 
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A. 3. The Dirichlet series J2ki>i {k^+py • Note that for our purposes, we only need 
g{n, 2a, 2b), so the series we are interested in is 

Z{a,h-s):=Y,Jj^r^fIy- (36) 
Clearly, this is closely related to the following function: 



Z^{a,b;s) := ^ 



(fc,Oez2\{o} ^ ^ 
= 4 ■ Z{a, b;s) + 2-[b = 0] C(2s - 2a) + 2 ■ [a = 0] ({2s - 2b) . (37) 

Informations on the poles and residues of Z{a, b; s) could be directly derived from 
the work of Pierrette Cassou-Nogues [21 P- 4 Iff], but there is a simpler way by using 
the reciprocity law for Jacobi's theta function. This reasoning is a generalization of 
Riemann's representation (see (HI section 1.12, (16)]) of ({s), 

TX-''^V{s/2) C{s) = + / + f /2) r'uj{t) dt, (38) 

s — 1 Ji 

where 

oo 

uj{t) = Y, e""'"* = 1/2 (^3(0, it) - 1) . (39) 

n=l 

Here, 6*3 denotes (one variant of) Jacobi's theta function (see j?! section 13.19, (8)]) 

00 

e,{z,t)= J2 ^"'e^™'' (40) 

n=— 00 

where q = e™^. 

A. 3.1. Jacobi's theta function. We rewrite ()4U|) by setting t = y and z = irx for x, y G C; 
i.e.: 

00 

^{x,y) = e'"K-"+f"') (41) 

n=— 00 

This series is absolute convergent for all x and all y with '^{y) > 0. Therefore, for fixed 
y the function f : z ^ •d{z,y) is an entire function. We have the properties 

i9(x + l,y) = 'd{x,y) , 

e2-i(-+f)^(x + y,?/) =^9(x,i/), 

which in fact determine the theta function up to a multiplicative factor c{y) (see |TH 
section 2.3]). Moreover, we have the following functional equation (see Theorem 
2.12, equation (2.29)]): 

^{x,y)=d(-,--\ e---4Jl 

\y yj \1 y 

Setting -di^y) := 19(0, iy) , we obtain as a special case the following reciprocity law, valid 
for all y with ^{y) > 0: 



(42) 
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^{y) is a holomorphic function in the half plane ^{y) > 0, with 3?(?/) = as essential 
singular line, see O Satz 2.13]. 

Interchanging summation, differentiation and integration in the appropriate places, 
we obtain 

-i [s) Z^:[a, b; s) = i [s 



{fc,/)ez2\{o} ^ ^ 



/•oo 

V / u'-' {-nky {-nl^Y e-^'^'+^'Hu (43) 

^^2\ rm "''0 



(fc,z)ez2\{o} 



1 / cl _^fc2„\ / a 



For (gni), we used T{s) = t'-^e-^dt = a' u^-^e-^^du with a = (P + ^2) 

Clearly, Z^,(a, 6;s) = Z^:{b,a;s). So w.l.o.g. we may assume a > b. We have to 
distinguish the following two cases, where we assume a > and b > 0: 

1 r°° 

—T{s)Z,{0,0;s)= / u'-^ {d{uf - 1) du, (44) 

t^r(s)Z.(a,M) = fu^-' {{£-A-)) {^A-))) (45) 

A. 3. 1.1. Case 1. Considering (j44j) . we basically repeat the reasoning in [TH p. 203]. 
From we get ^'*9(^)^ = 'd{uf' . Assuming 3?(s) > 1, we compute: 



00 



n^^^ (?9(u)^ - 1) du 



1 1 POO poo 

= — + -+ / {^{tf - 1) dt + / n^-i (^^(u)' - 1) dn. (46) 

^ s — 1 Ji Ji 

The integrals in the last line converge for all s G C and constitute holomorphic functions, 
so Z^{0, 0; s) only has a simple pole for s = 1. 

A. 3. 1.2. Case 2. Considering we have to adjust the preceding method appropri- 
ately. For convenience, set da{u) := ■^'&{u) . Then we have: 

/» 1 roo 

u-'~' (Mu) M^)) du = j u'-' {Mu) M^)) du + j u-'-' (4(n) M^)) du. 

Now use fl42j) in the form 
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and combine this with the the formula 
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y \y 



(2a)- 



A:=0 



la—k 



dy 



a—k 



f 



1 ^ yk~2a~^ 

y ■ 



to obtain 



-1 



^,+,^^(2a)^(26) 



EE 

fc=0 j=0 



4'^+''a!6! 

2fc ^r,.N2i 



4^+3 k\j\ 



U \U 



[a = k M) = j] \ u 



k+j-2a-~2b~l 



Now in the same way as before, this gives 



(-1)"+'' (2a)!(26)! 
4'^+''a!fe! {s-a-b-1 

(2a)- (26)^ 



EE 

k=0 j=0 



4k+jk\j\ 



{^a-k{u) ^b-j{u) -[a = kAb = j]) (47) 



Again, the integrals converge for all 2; G C and constitute holomorphic functions, whence 
Z^{a, b; s) has only one simple pole ats = a + 6+ l. 

We summarize all this information in the following proposition. 

Proposition 3. For arbitrary nonnegative integers a, b, the series ;>i (^k^'^py is 
convergent in the half-plane 9ft(z) > a + b + 1 and defines a meromorphic function 
Z{a, b; z) : C ^ C with a simple pole at z = a + b+1, where the residue is ia-yb+i^lfl'a+by. ■ 
If a > and b > 0, this is the only pole. 

If a = (or b = 0), there is another simple pole at z = b + \ (or z = a + \), where 
the residue is —\ ([6 = 0] + [a = 0]). 

Moreover, we have the following information on special evaluations of Z{a, b; z): 



Z{a, b; —n) = -[a = 6 = n = 0] for n = 0, 1, 2, 



(48) 
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For the absolute term Ca^b in the Laurent series expansion of Z{a, b; z) at zq = a + b+ ^, 
we have the following formulas: 

1 - 
Co,o = -7-l + 2 J t-^^{tf-l)dt, 

7 1 A''-^a\((-lT + 1) r°° 1- 

c . - ^"^'"(" + ^)' /^_2M1M! + r_ir+^ t^+^-i^ (t) A(t) dt 

(2a + 26)! 1^ 4^+^a\b\ + ^ A ^ ^^a(t) dt 

A^(2a)^W 
fc=0 j=o '-^ ' 

X / (^^a-fc(t) ^b^j{t) -[a = kAb = j]) t'^+^-'^^^-sdt ) /or a > 6 > 0. (49) 



Proof. This information is extracted straightforwardly from (p?7j) together with ()44j) . 
(jini) and (gni), (HZI), respectively. 

(The evaluation (^(— 2n) = for nonnegative integers n, which is needed for fl48p . can 
be found in 1.13, (22)].) □ 
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